Spin drag in noncondensed Bose gases 
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We show how time-dependent magnetic fields lead to spin motive forces and spin drag in a spinor 
Bose gas. We propose to observe these effects in a toroidal trap and analyze this particular proposal 
in some detail. In the linear-response regime we define a transport coefficient that is analogous 
to the usual drag resistivity in electron bilayer systems. Due to Bose enhancement of atom-atom 
scattering, this coefficient strongly increases as temperature is lowered. We also investigate the 
effects of heating. 
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Introduction: Coulomb drag and spin Coulomb drag — 
Understanding electronic transport [Ij is one of the most 
important goals of condensed-matter physics. Indeed, 
materials are often characterized according to their trans- 
port properties. Furthermore, transport measurements 
provide important physical information. For example, 
the temperature dependence of transport coefficients, like 
resistivity and conductivity, contain information on the 
elementary excitations and their scattering mechanisms. 
Moreover, the magnetic-field dependence allows for ex- 
tracting the electronic phase-coherence length. 

Analyzing results of transport measurements is com- 
plicated by the multitude of effects, like electron-electron 
and electron-phonon interactions, that contribute. This 
problem is to a large extent circumvented in the Coulomb 
drag measurement of Gramila et al. 0, illustrated in 
Fig. [U that aims at singling out the electron-electron 
interactions from the start. In this setup a bilayer of 
two-dimensional electron gases is separated by a tunnel 
barrier. A current / is driven through the bottom layer 
that drags along the electrons in the other layer. In the 
top layer, an electrochemical potential is built up that 
cancels the drag and induces a voltage drop Vd, which 
results in a drag resistivity pn = Vd/I- Originating from 
electron-electron interactions, this resistivity usually has 
the typical pjj oc T 2 Fermi-liquid- like temperature de- 
pendence at low temperatures 0, 0, @] . 

One approach in describing Coulomb drag is in terms 
of a function Trj (v t — Vb) that gives the rate of change 
of momentum per volume of the electron gases due to 
Coulomb scattering 0, H] , and that to a good approxi- 
mation depends on the difference in the drift velocities 
v t and Vb in the top and bottom layers. The equations of 
motion for these drift velocities are then given by 



FIG. 1: Illustration of a Coulomb drag measurement: a pair 
of two-dimensional electron gases is separated by a tunnel 
barrier. A current is applied in one of the layers leading to 
a voltage drop in the other. (Adapted with permission from 
http://www.sp. phy.cam.ac.uk/SPWeb/| ) 



in the top (bottom) layer are denoted by nt(nb) and 
E t (Eb), respectively. We have also added a scattering 
rate 1/r* and 1/tj, for the top and bottom layer, that ef- 
fectively takes into account intralayer Coulomb scatter- 
ing, electron-phonon interactions, and disorder. In ap- 
plying the above result to the situation in Fig. [T]we take 
v t = since there is no current in the top layer. Solving 
for Vb in the steady state, using that <C nbm e Vb/Tb 
because the interlayer scattering is much weaker than the 
intralayer scattering, we find that Vb = —eEbTb/m e , as 
usual. In the linear response regime the drift velocities 
are small, and we can use that Tu(v) ~ T' D (0)v incorpo- 
rating the fact that there is no net momentum transfer if 
the drift velocities are equal. We then find that the elec- 
tric field in the top layer is E t = — T' D (0)vb/en t . Using 
that the current density in the bottom layer jb = —ribevb, 
we have for the drag resistivity 



n t m e 



dv t 
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-T D (v t - v b ) - n t eE t 
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PD = ~ 
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(2) 



nbm e Vb 

nbm e —— = +TD(vt - Vb) - meEb . (1) 

at Tb 

Here, — e is the charge of the electron and m e its mass. 
Furthermore, the electronic density and the electric field 



This result shows that the drag resistivity is determined 
by the slope of the function Trj(v) at v = 0. 

In an analogy to Coulomb drag, D'Amico and Vignale 
proposed spin Coulomb drag 0, which was observed by 
Weber et al. Q ■ Spin drag, in which the layer degree of 



freedom from Coulomb drag is played by the spin of the 
electrons, is very similar to Coulomb drag. In this Letter 
we study spin drag due to the short-range interatomic 
interactions in a spin-one Bose gas in the normal state, 
and propose an experiment to observe it making use of 
so-called spin motive forces. For this system, we derive 
equations of motion similar to Eqs. (Q}. The absence of 
disorder and an underlying lattice that supports phonons 
implies that the analogues of the scattering times t& , r t 
are infinite. Nonetheless, we recover a great deal of the 
phenomenology of conventional electronic transport. In 
particular, we define a transport coefficient analogous to 
Pd which for bosons becomes large at small temperatures 
due to Bose enhancement, i.e., the enhanced scattering of 
bosons to states that are already occupied. In addition, 
we investigate heating effects and find that they are com- 
pletely analogous to the usual Joule heating in electronic 
systems. 

Ultracold atomic gases and spin motive forces — We 
consider ultracold atoms with hyperfine spin F in a time 
and position dependent magnetic field with a direction 
given by the unit vector fi(x, t), such that the Zeeman 
interaction reads Hz = — Afl(x, i) ■ F/h, where F are 
the spin operators and A is an effective Zeeman splitting 
energy. If the magnetic-field direction is varying slowly 
in space and time, it is convenient to choose S~2(x, t) as 
the local spin quantization axis. In this frame of refer- 
ence spatial and temporal variation of the magnetic-field 
direction manifests itself as fictitious, or fixed-frame, elec- 
tric and magnetic fields E and B [8[ that are ultimately 
due to the spin Berry phase Q. For atoms with spin 
projection mp these arc given by 

'on 



B„ 



= mpMl 



dt 



x V q O 



mphe a /3 7 £l ■ (V/3O x V 7 f2) 



(3) 



where h is Planck's constant, e Q/ 3 7 is the three- 
dimensional fully antisymmetric Levi-Civita tensor and 
a sum over repeated Cartesian indices a, (3, 7 6 {x, y, z} 
is implied. Note that, because the atoms are neutral, 
there are no real electromagnetic fields that couple to the 
atomic motion. In the context of ferromagnetic metals 
these fictitious electric and magnetic fields respectively 
underlie the phenomena of spin motive forces induced 
by moving domain walls and the topological Hall effect, 
both of which have been observed very recently 1(J 11 1 • 
In the context of cold atoms, the Aharonov-Bohm phase 
due to the fictitious magnetic field, in combination with 
phase coherence, has been used to imprint coreless vor- 
tices on F — 1 spinor Bose-Einstein condensates 12, 13j |. 
For the existence of the fictitious electric and magnetic 
fields phase coherence is, however, not required [8[ and 
we can focus instead on the semi-classical regime using 
the equation of motion 

m dv ™ F = F mp + v mp x B mp , (4) 
at 




FIG. 2: (Color online) Illustration of spin drag in a toroidal 
trap. The atoms with ttif = 1 (red circle with arrow) are 
accelerated by the motive force due to the time-dependent 
quadrupole field (blue arrows) of the IofTe-Pritchard trap. 
Due to interactions the atoms with mp = are dragged along. 



for an atom with velocity v mF and spin projection mp. 

The specific geometry we consider is illustrated in 
Fig. [2] and consists of a toroidal trap with radius R and 
effective cross section area A in the transverse direction, 
created by a rapidly- moving laser beam To im- 

plement the fictitious electric field we superpose a Ioffc- 
Pritchard magnetic trap. Fictitious electric fields along 
the torus are achieved by varying the axial bias field of 
the Ioffc-Pritchard trap, so that [13[ 



si{<t>) = n z (t)z + y/i-w z (t) 



f cos 2(f> — sin 20 



(5) 



in cylindrical coordinates (r, 
find that 



z) . Using Eq. © we then 



E,, 



2hmi 



1 dft z (t) 
R dt 



(6) 



and B, 



0. The adiabatic approximation that leads 



to the above holds when the timescale To on which the 
direction of the external magnetic field is changed is much 
larger than the spin precession time h/A. Furthermore, 
this spin precession time should be smaller than the time 
it takes the atoms to encircle the torus. Since A is a large 
energy scale, these conditions are easily satisfied. 

Spin drag — We now specifically consider noncon- 
densed bosonic atoms with F = 1, e.g., sodium atoms. 
Furthermore, to study spin drag we consider the case that 
the trap is loaded with equal densities of atoms in spin 
state |1) with mp = +1 and |0) with mp = 0. According 
to Eq. (JHJ) the atoms in spin state |1) then feel a fictitious 
electric field E± = E along the torus which accelerates 
them. The atoms in the |0) state feel no fictitious elec- 
tric field but may accelerate due to spin drag, i.e., due to 
collisions with the other atoms. 

To investigate the spin drag quantitatively, we use an 
effective one-dimensional Boltzmann equation for the dis- 
tribution function fi(k, t) and fo(k 7 1) of the mp = 1 and 
mp = atoms, with hk the momentum along the torus, 
given by 



3 



, Ed ft _ 2tt , 2B ,2 f dk 2 f dk 3 f dk 4 

x{[l + / 1 (fc,i)][l + / (fc 2 ,i)]/ 1 (fc 3 ,t)/o(fc4,t)-/l(fc,<)/o(fc2,t)[l + A(fc 3 ,t)][l + /o(fc4,t)]} , (7) 



and where the equation of motion for fo(k 7 t) follows by 
taking E = and interchanging fo(k, t) and fi(k,t). The 
interspecies collisions are determined by the two-body T- 
matrix — 4irah 2 /mA, with a the scattering length 
for collisions of atoms between atoms in states |1) and |0), 
and m the atomic mass. The single-particle dispersion is 
£k = h 2 k 2 /2m. On the right-hand side we have ignored 
intra-spin-species collisions which tend to restore local 
equilibrium and are zero in the approximations outlined 
below. Also note that, contrary to electronic transport in 
solid-state physics, there are no terms corresponding to 
elastic or electron-phonon collisions because in cold-atom 
systems there is no extrinsic disorder or an underlying 
ionic lattice. 

Since the intraspccies collision enforce local equi- 
librium for each spin species we use a Bose-Einstein 
distribution function with nonzero drift velocity as 
an ansatz to solve the above equation. Specifically, 
we take fi(k,t) = N B (e k _ mvl{t)/h ) and /o(M) = 
N B (e k - m v {t)/h), with7V B (e) = [ e /M*-/0_i]-i the Bosc- 
Einstein distribution function at chemical potential fi and 
inverse thermal energy (3t = 1/fcsT. In first instance we 
take the temperature constant in time. The time de- 
pendence of the chemical potential is determined by the 
conservation of the number of atoms in each spin state 
and is left implicit. From the Boltzmann equation we 
then find that, cf. Eqs. ([1]), 

nm—^- = nE + T(vq — v i) ; 
at 

nm—— = -r(u -i>i), (8) 
at 

where n is the one-dimensional density, and the func- 
tion that determines the rate of momentum transfer from 
species |1) to |0) is found from Eq. ([7]) as 

/dk 
—T coll (k)hk , (9) 

with the right-hand side evaluated using the shifted Bosc- 
Einstcin distribution functions. In Fig. [3] we plot this 
function for various values of the degeneracy parameter 
nA, with A = y27r ft 2 /mfeT the deBroglie wave length. 
We find that in the classical limit nA — > it is given 
by T(v) = (4T:ahn) 2 Erf(mAv/h)/A 2 m. For increasing 
degeneracy T(v) develops local maxima and minima at 
small |w max | which are due to Bose enhancement of inter- 
species scattering. 

From the equations of motion in Eq. we see that the 
sum of drift velocities increases indefinitely. The relative 
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FIG. 3: Rate of momentum transfer F in units of 
A 2 m/ (2ahn) 2 . The difference of drift velocities is in units 
of h/mA. The inset shows the derivative of F at v = 0. 



drift velocity v = v\ — vq can approach a steady state, 
provided the motive force E is not too large. That is, 
from Eq. we find that if nE < 2r(w max ) the sys- 
tem approaches a steady state where dv/dt = 0. In 
the linear-response regime E and v are small and we 
have that T(vq — Vi) ~ r'(0)(vo — ^i). Introducing the 
relative-momentum particle current j = n(vi — vq), we 
have in linear response that Vi — vq = nE/2T'(0). From 
this we define in the linear-response regime a resistivity 
p = E/j = 2r'(0)/n 2 , that is analogous to the drag re- 
sistivity in Eq. For fermionic atoms this resistivity 
would vanish at small temperatures. For bosons it be- 
comes larger due to Bose enhancement. This is further 
illustrated in the inset of Fig. [3] which shows r'(0) as a 
function of l/(?iA) 2 cx T. 

The fact that the total kinetic energy of the system is 
increasing suggests that beyond-linear-response effects, 
such as heating, may be important. To investigate these, 
we have to solve Eqs. ([8]) coupled to an equation for the 
temperature. This equation is most easily derived by 
considering the total energy U = J(fi + fo)ekdk/2TT. We 
evaluate this energy within our ansatz of Bose-Einstein 
distribution functions with nonzero drift velocities and 
time-dependent temperature T(t) in this case. Using 
the Boltzmann equation in Eq. and Eqs. (JSJ), we 
find that dQ/dt = n[v\ — Vq)T{v\ — vq), where Q 
I ' — mn(v 2 + t>o)/2. This energy is determined by the 
spread in velocities in the gas of atoms and is therefore a 
measure for its temperature. Evaluating the above using 
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FIG. 4: Drift velocities of as a function of time, for nA = 45 
(solid lines) and nA = 30 (dashed lines). The upper and lower 
curve correspond to «i and Vo, respectively. 

the linear-response expression for the difference in drift 
velocities, we find that dQ/dt = npj 2 /2, which is analo- 
gous to Ohmic heating in electronic systems. 

We go beyond linear response by solving the equation 
for dQ/dt coupled to Eqs. ©. We consider the case 
that the axial magnetic field of the Ioffe-Pritchard trap is 
inverted in a time T , so that Q z (t) = (2t— 1)/T for < 
t < To, and constant for t > Tq. This implies via Eq. (J6]) 
that E mF = Ah4>/RT for < t < T and zero for t > T . 
We consider specifically 23 Na atoms. As parameters we 
take T = 10 ms [H, R = 5 /m, T = 400 nK. For the 
one-dimensional density we take n = 10 12 cm -3 x A, with 
A = 7r (5^m) 2 [15| . For these parameters nA = 45. The 
result is shown in Fig.rj] together with the result for nA = 
30. We find that heating effects are negligible on the time 
scale that is shown. For each pair of curves the upper one 
corresponds to v\(t), which, due to acceleration by the 
motive force, acquires the value v\ ~ Ah/mR in the time 
To. The lower curve corresponds to Vq which starts at 
vq(0) = 0. Due to the spin drag, the latter velocity also 
becomes nonzero, which can be experimentally measured 
by studying the momentum distribution after expansion. 
Note that the drag effect is larger for larger nA due to 
the Bosc enhancement. 

Discussion and conclusions — There are other exper- 
imental setups possible to observe spin drag effects. For 
example, a cigar-shaped optical trap together with a 
magnetic field gradient in the axial direction also leads to 
relative motion of the two spin species. We note that the 
hydrodynamic regime, where spin-drag effects should be 
large, has been realized recently in such cigar-shaped sys- 
tems (l6j . Another possibility is using a Raman transi- 
tion to convert a fraction of the atoms of one spin species 
to another, and to set them into motion with a velocity 
determined by the recoil energy of the two-photon pro- 
cess. Such an experiment has already been performed 
with Bose-Einstein condensates [13] ■ However, to study 



spin drag, and other analogues of electronic transport, 
the noncondensed case is more suitable because the in- 
coherent collisions of the thermal atoms, rather than the 
coherent interactions in a Bose-Einstein condensate, are 
analogous to the collisions of the electrons. 

Yet another experimental possibility is to use a 
sinusoidally-varying axial bias field of the Ioffe-Pritchard 
trap. This would lead to an AC electric field, and the pos- 
sibility to measure the frequency dependence of the trans- 
port coefficient p. Other interesting generalizations of the 
present work are including mcsoscopic phase-coherence 
effects, and effects of critical fluctuations. Drag effects 
can also be measured in Fermi gases, leading to another 
way to probe the many-body physics of these systems. 
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